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TECHNICAL NOTE U067 


APPROXIMATE ANALYSIS OF EFFECTS OF LARGE DEFLECTIONS 
AND INITIAL TWIST ON TORSIONAL STIFFNESS OF A 
CANTILEVER PLATE SUBJECTED TO 
THERMAL STRESSES 

By Richard R. Heldenfels and Louis F. Vosteen 

SUMMARY 


An approximate analysis of the no nlin ear effects of Initial twist 
and large deflections on the torsional stiffness of a cantilever plate 
subjected to a nonuniform, temperature distribution iB presented. The 
Von K£rm£n large-deflection equations are satisfied through the use of 
a variational principle. The results show that initial twist and applied 
moments can have significant effects on the changes in stiffness produced 
by nonuniform heating, particularly in the region of the buckling temper- 
ature difference. Results calculated by this approximate analysis are 
in satisfactory agreement with measured torsional deformations and changes 
in natural frequency. 


INTRODUCTION 


One of the structural problems of high-speed flight is the reduction 
of effective stiffness of structures due to the thermal stresses produced 
by aerodynamic heating. A reduction in torsional stiffness can be an 
important factor in aeroelastic problems as indicated in references 1 
and 2. A similar reduction in stiffness produced by thermal, stresses 
presumably caused the flutter and failures of seme structural models 
described in reference 3 . A simple method for calculating the reduction 
in torsional stiffness of thin wings is presented in reference 4. .In 
reference 5 the results calculated from a small-deflection plate theory 
are compared with experimentally determined changes in the torsional 
stiffness of a cantilever plate rapidly heated along the longitudinal 
edges. The theory used in reference 5 predicted the general effect of . 
thermal stresses on the torsional stiffness, as indicated by measurements 
of torsional deformation and changes in natural frequency of vibration, 
but overestimated the magnitude of the changes. 
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The purpose of this paper is to present the results of an approximate 
analysis to show that the differences between theory and experiment noted 
in reference 5 are due to the nonlinear effects of large deflections and 
initial deformations not included in the small-deflection analysis. The 
analytical approach used to account for large deflections and initial 
deformations is presented, the general significance of the results is 
discussed, and calculated values are compared with the experimental data 
of reference 5 . 


SYMBOLS 


• • • A-j_ 

b 

C 

D 

E 

F 

fl, f 2 , . . .f-j_ 

I l> 1 2 ' • * •■’•i 
GJ 

(GJ ) 0 

A(GJ) 

aCgj^ 

aCgjW 

m 

M 


plate length in x-direction 

coefficients of series expansion for plate deflection 


half -plate width in y-direction 
coefficient of stress function 


plate flexural stiffness. 



) 


modulus of elasticity 

stress function defining stress distribution in plate 
selected functions of x and y 


values obtained from definite integrals 


torsional stiffness 

torsional stiffness of flat, unstressed plate, 
incremental torsional stiffness 

initial incremental torsional stiffness before heating 
minimum incremental torsional stiffness 


nondimensional moment 

applied moment (positive in direction of positive twist) 
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T 

AT 


AT 


cr 


w 


a 

e 

i 

0 

e i 

A 

^ min 

P 

P 

cr x-»°y 


T xy 


T 


<P 


'Pi 

cu 


pressure 

plate thickness 

temperature 

particular temperature difference 
critical value of AT 

total plate deflection 
initial plate deflection 
coordinate axes 

coefficient of thermal expansion 
small dynamic perturbation 

exponential parameter of temperature distribution function 

twist at plate tip 

initial twist at plate tip 

temperature ratio, AT/AT cr 

value of A corresponding to A(GJ) m -j n 

Poisson's ratio 

density 

normal stresses in plane of plate in x- and y-directions , 
respectively, positive for tension 

shear stress in plane of plate 
time 

nondimensional twist 
initial nondimensional twist 
circular frequency 

circular frequency of unheated perfect plate 
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ay. 

«kLn 

V 2 


initial circular frequency 
minimum circular frequency 

differential operator, + 2_ 

dx 2 by £ 


differential operator. 


+ 2 


Sx 2 dy 2 Sy^ 


Double dot indicates second derivative with respect to time. 


ANALYSIS 


Statement of Problem 


The studies presented herein are primarily concerned with the twist 
of a uniformly thick rectangular cantilever plate shown in figure 1. The 
equations are derived to consider the effects of initial twist, applied 
moments, and thermal stresses on the torsional deformations and natural 
frequencies of plates. The analysis involves an approximate solution of 
the Von K&rmAn large-deflection equations which have been modified to 
include the effects of initial imperfections and nonuniform temperature 
distributions. The modified equations are (from ref. 6) 


4 

V F = 


-Ea 




+ E 


7 c) 2 w \ 2 

d 2 w 3 2 w 

^ 2 w t \ 2 a 2 *! 

S 2 w ± 

\Sx By/ 

dx 2 3y 2 

W by j dx 2 

dy 2 


(la) 


(w - w-i ) =6 

, t fd 2 F b% . b^F b% 

0 

•a 2 w\ 

(lb) 


u 

\Sy 2 Sx 2 bx by 

dx by 

bx by j 

where F is 

the stress 

function such that 
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(2a) 
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°y 


S 2 f 

dx^ 


(2b) 


T 


xy 


b 2 ? 

Sx by 


(2c) 


Equation (la) is the differential equation for compatibility in the 
plane of the plate and equation (lb) is the differential equation for 
equilibrium of forces acting perpendicular to the plane of the plate. 

The solution of equations (l) , subject to the proper boundary conditions 
on the deflections and stresses, describes the behavior of the plate under 
the applied loads p and the temperature distribution T. For the pur- 
pose of this paper the assumption is made that no external loads are 
applied in the plane of the plate so that the inplane stresses defined 
by the stress function F are zero at the free boundaries (stress-free 
edges) . Equations (l) apply to both the dynamic and the static problem 
if p includes the dynamic loads as well as static loads. 


Approximate Solution of Equations 


Exact solutions of the Von KArmAn large-deflection equations are 
difficult to obtain, but approximate solutions can be obtained to any 
desired accuracy by several methods. For the analysis herein equations (l) 
will be satisfied through the use of the following variational principle: 




(5) 
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This principle is essentially a particularization of one given hy 
Reissner (ref. 7) and also has teen modified to include the effects of 
initial deformations and a nonuniform temperature distribution. The 
deflection w must satisfy the geometrical boundary conditions, that 
is, the conditions imposed on slopes and deflections. Variation of 
equation ( 3 ) with respect to F yields equation (la) and the associated 
natural boundary conditions, and variation with respect to w yields 
equation (lb) and the associated natural boundary conditions. 

The following assumptions for stress, deflection, and temperature 
are made : 


F = Cf x 
w = b9fg 
w^ = b0j_f2 
T = ATf 5 


(4) 


where the stress coefficient C and the tip twist 8 are to be deter- 
mined by means of the variational principle. The functions f j_, f 2 , 

and f^, the initial tip twist 0^, and the temperature coefficient AT 
are presumed to be known. Selection of these quantities is discussed in 
appendix A. It should be noted here that, although the temperature 
distribution is considered to vary with time, the shapes of the plate 
deflection and of the stress function are assumed to remain constant 
during heating or loading. 

The unknown coefficient of the stress function C is obtained by 
substituting equations (4) into equation ( 3 ) and taking the variation 
with respect to C. The result is 



where 



- 2p- 


a 2 fi a 2 fi 

dx 2 By 2 


+ 2(1 + 11) 



dx dy 
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The relation between the twist, the load, and the temperature is 
obtained by performing the variation on equation ( 3 ) with respect to the 
twist 0 with C held constant. The resulting relationship is 



®±)ij+ 


2tcb 2 ei, 


- b iT 


pf 2 dx dy = 0 


(9) 


where 







O Io o JL o , x 

2p— ^ + 2(1 + p) 

ax 2 dy 2 



dx dy 


( 10 ) 


The pressure loading p now will be considered to consist of a static 
loading p s and a dynamic or inertia loading pt so that 


P = P s - Pt w 


(H) 


For this analysis the lateral static load p s will be restricted 
to two equal concentrated loads P applied at the comers x = a, 
y = ±b in ‘such a way as to form a couple about the x-axis. Then 


1 


Pf 2 a* dy = JLi 5 . ptb'eig 


( 12 ) 
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where 


I5 = JJ S(x-a) [s(y-b) - 5 fg dx dy 


*-JT 


fo dx dy 


(13) 

(14) 


M = ZbP - (15) 


and 8 is the Dirac delta function. 

If the value of C from equation (5) is substituted into equa- 
tion (9) and the terms are rearranged, the resulting expression is 


l > 2 ( 9 - 0 i) - §15 


tb 

D 


0I 6 + 


2tb 2 8 


Eh, AT 


I3I2 

In 


Eefe 2 

d ' 



(16) 


which gives the relationship "between twist, applied moment, inertia 
loading, and temperature . 


Definition of Parameters 

It is convenient to define certain parameters which can be obtained 
from equation (l6) when the large-deflection effects are neglected and 
the initial twist 0^ is 0. Then, from equation (l6) 

b 2 9!4 = JJjI 5 - p ^ 0I 6 + Ex AT Sfe (17) 


If the plate is not vibrating and no loads are applied, the critical, 
or buckling, temperature difference is obtained from equation (17) as 

AT..,, , -S- SI 
or 2Bxt I 2 I 3 


(18) 
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Similarly, when the plate is not vibrating and the temperature is 
uniform (AT = 0) hut the load is acting, the moment-twist relation is 
given by 


a _ Ma 
(GJ ) 0 


( 19 ) 


where (GJ)q is the torsional stiffness of the flat plate given by 


(GJ) 0 = 2ab 2 D 


( 20 ) 


In the absence of heating and loading the frequency may be found to be 


<*0 


2 _ D_ 

P* x 6 


( 21 ) 


When equations (l8), ( 19 ) , (20), and (21) are substituted into equa- 
tion ( 16 ) the result may be written as 


e.ej.JV 1 S + ^e-s( 9 2 - 01 2 ) 
(gj) 0 ' ' 


2\ 2Etb c 


D 


M! 

Vk 


(22) 


If a "nondimensional" twist is defined so that 


9 = 


2Etb 2 fe) 2 
D I X I 4 


(23) 


along with a nondimensional moment 


Ma 2Etb 2 j 1 ?) 2 

(oj) 0 y d i x i4 


m 


and a temperature ratio 
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equation (22) may be written simply as 


cp - cp^ = m + Acp - cpfcp 




(26) 


Application to Static Problem 

The relationship between twist, moment, and temperature is given by 

equation ( 26 ) when the dynamic term cp - is 0 . Equation (26) then 

0^2 

becomes 


cp-q^sm + Acp-cp 



( 27 ) 


The incremental stiffness of the plate, defined as the rate of change of 
moment with respect to twist, is given by 

A(GJ) _ a 3 m 

(gj ) 0 (gj ) 0 de 


or, from equation (27), 

(£8) 

Application to Dynamic Problem 

In order to determine the effect of temperature, moment, and Initial 
twist on the natural frequency of torsional vibration, the quantity cp 
in equation (26) is replaced by 9 e sin cjcfT where 9 is considered 
to be the static solution obtained from equation (27) and e represents 
a small dynamic perturbation about the static equilibrium position. Sub- 
tracting equation (27) from the perturbate'd relation, neglecting higher 
order terms in e, and dividing by e Bin orr yields 
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. The frequency may then he written in terms of the twist and tempera- 
ture as 

= 1 - A + 3 q£ - fRj 2 (29) 

which is identical to the incremental stiffness (eq. (28)). This identity 
results from assuming similar deflection mode shapes for thermal buckling, 
twist due to applied moment, and torsional vibration and would not be 
expected to apply if the mode shapes involved were different. 


RESULTS AND DISCUSSION 


General Results of Equations 

The behavior of the plate as determined by equation (26) is discussed 
in the following sections for several combinations of conditions . In most 
cases the calculations include large values of initial twist and applied 
moment and have been carried well into the region where the temperature 
difference AT exceeds the critical value. Although some of the results 
may be beyond the range for which the analysis is accurate, these results 
have been presented to illustrate the trends indicate^ by the equations 
despite the fact that they may not be quantitatively correct in some 
regions . 

Twisting due to an applied moment . - The twisting of the plate due 
to an applied moment is given by equation (27} if the temperature ratio A 
is set equal to 0 . The results from this expression are presented fig- 
ure 2 for various values of the initial twist. 

In figure 2 , small-deflection results would plot as lines at 45 ° to 
the coordinate axes. The large-deflection results become increasingly 
different as the moment or initial twist is increased and indicate that 
the plate becomes substantially stiffer as the twist is increased. 

Buck! lug due to nonuniform heating . - The twist of the plate (as sinning 
that the buckling mode is a twisting action) is given by equation (27) 
when the moment m is 0; for the initially flat plate cp. would also 
be 0 . The initially flat, or perfect, plate begins to deform only after 
the critical temperature difference is reached, whereas, as indicated by 
equation (27) , the plate with initial twist begins to deform immediately 
upon heating. The results obtained from evaluation of equation (27) are 
given in figure J for several values of the initial twist. 
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If the initial twist is large, a plot of twist against temperature 
ratio, like figure 3, does not give an accurate indication of the huckling 
temperature . Only if the initial twist is small is there a definite knee 
in the curve as the huc kli ng temperature is approached, hut this knee 
occurs "below the huckling temperature of the perfect plate. 

Combined action of applied moment and nonuniform heating .- Equa- 
tion (27) applies directly to this case hut may he more conveniently 
written as 


cp - cp> 


A'cp + cplcp 



(30) 


where 


cp' = qijL + m 


(31a) 


and . 


A' = A 




(31h) 


Equation (30) now has the same form an equation (27) when m is set 
equal to 0. The results plotted in figure 3 then apply to this case 
also if cp^ is replaced "by cp' and A hy A' . The moment then effec- 
tively acts the same as an initial twist if" the temperature ratio is 
reduced hy the quantity m 2 + 2mcpj_. 

In figure 4 the relation between moment and twiBt has heen indicated 
for two values of initial twist and several values of A. These curves 
show the characteristic increase in stiffness as the twist increases hut 
also show a reduction in stiffness as the temperature ratio increases. 
These changes in stiffness are examined further in the following section. 
The portions 'of the curves where ;the slope is negative have heen shown 
as dashed lines and are regions of unstable equilibrium which would not 
exist in the physical problem. The unstable portion of the curve exists 

whenever A + cp^ 2 - 3q£ > 1. 

Frequency and incremental stiffness . - As has heen noted, the square 
of the frequency ratio (eq. (29) ) varies in the same manner as the incre- 
mental stiffness (eq. (28)) and, therefore, any of the following discus- 
sion pertaining to stiffness applies directly to the square of the fre- 
quency. Also, the figures which are presented for the incremental 
stiffness have heen labeled with the square of the frequency ratio as 
well as with the stiffness ratio. 
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The variation in incremental stiffness with temperature given "by- 
equation ( 28 ) is shown in figure 5 for several values of initial twist 
when m is equal to 0. The results show that, when the perfect plate 
is heated, the stiffness decreases linearly as the thermal stresses 
develop and becomes zero when the buckling temperature difference is 
reached. Further heating causes the plate to twist and the twisting in 
turn leads to an increase in stiffness at a rate twice that of the initial 
decrease. 

If the plate has an initial twist, the heating causes the twist to 
increase as the heating progresses. This twist leads to an increase in 
stiffness that tends to counteract the reduction produced by thermal 
stresses and, as a result, the stiffness first decreases and then 
increases without going to zero. The points of minimum stiffness occur 
at temperatures lower than the buckling temperature. The locus of points 
of minimum stiffness has been indicated in figure 5 hy the dashed line. 

If the initial twist is sufficiently large, the stiffness does not 
decrease but increases with heating. 

When a moment is applied to the plate, figure 5 will indicate the 
incremental stiffness of the plate if qx is replaced by cp 1 and A 
by A 1 (as defined by eqs. (31a) and (31b)). 

The locus of points where the incremental stiffness Is a minimum 
is given by the equations 



Equations ( 33 ) and (3*0 show that the nonuniform heating will increase 
the stiffness of the plate if the initial twist is greater than that 
given by 




14 


NACA TN 4067 


The initial incremental stiffness of the plate is given hy equa- 
tion (28) when A is 0. The ratio of the minimum to the initial incre- 
mental stiffness therefore may "be written as 


A(GJ)i 1 + 3q^ - cpj. 2 


(35) 


Some results from equation (35) have "been plotted in figure 6. This fig- 
ure indicates that small changes in the initial twist may cause large 
changes in the minimum incremental stiffness. The effects of applied, 
moment are again similar to those of initial twist and certain combina- 
tions of the two can lead to. drastic stiffness changes. For negative 
values of cp^ the curves would he similar, with the stiffness ratio 
going to zero whenever cp^ = -m. 


Comparison With Experiment 

Results calculated with the previously derived equations are compared 
in the following sections with the experimental results reported in refer- 
ence 5- The results given in reference 5 are for a square cantilever 
plate which was heated along the two longitudinal edges hy carbon-rod 
radiators. Typical temperature histories of ^points on the heated edge 
and on the longitudinal center line are given in figure 7* Heat was 
supplied to the plate edges for 16 seconds; then the plate was allowed 
to cool. The deformations of the plate under the influence of nonuniform 
heating were determined for the conditions of no load and applied positive 
and negative tip moments. The changes in natural frequency of the first 
torsion mode were also measured during a heating test. 


The expressions used in reference 5 for stresses, deflections, and 
temperature are retained except that herein only three terms are used 
for the deflections. In nondimensional form these expressions are 


f l = 







+ A- 


■(f) + Mf) 


(36) 


1 + +• A 2 


(37) 
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and 


*3 . {if (38) 

The function f-j_ satisfies the condition of zero stress on the free 
edges and f 2 specifies zero slope and deflection along the root. Hie 

two undeter min ed coefficients of equation ( 37 ) are established from the 
small-deflection buckling analysis and carried through the large- 
deflection analysis as constants. Their values are given in appendix A 
along with data on the influence of the number of terms in the deflection 
function on the accuracy with which the buckling temperature difference, 
natural frequency, and ratio of twist to an applied moment can be calcu- 
lated. The exponential parameter £ of the temperature function is 
adjusted to appro xima te the variation in the temperature distribution 
during a test. 

The numerical evaluation of the integrals and related functions 
required for comparison of theory and experiment are given in appendix B. 

Initial plate shape .- The initial shape of the plate of reference 5 
was measured and is indicated in figure 8. In the selection of a value 
of initial twist, the higher order shapes (which are unlikely to have 
much influence on the twisting) have been Ignored. The free corners 
were connected by a straight line in order to obtain a value of Initial 
twist of 0 i = 0 . 35 ° which was used for comparison of theory and 
experiment . 

Twist due to an applied moment .- The deformations resulting from 
heating for three values of applied tip moment are presented in figure 9 
and compared with curves calculated by use of a value of 9-j_ = O. 35 0 . 

The agreement between theory and experiment is satisfactory, although 
the theory overestimates the twist in the vicinity of the maximum temper- 
ature difference (about 1 6 seconds) . No theoretical results are presented 
past 20 seconds because, beyond this time, the actual temperature distri- 
bution cannot be represented very well by the one-parameter temperature 
function. 

Another comparison of measured and calculated deformations is shown 
in figure 10 where the abscissa is the temperature difference instead of 
time. Because the shape of the temperature distribution changes with 
time, the portion of the curve for decreasing AT does not retrace the 
heating portion. 

Natural frequency .- The changes in natural frequency during a heating 
test are shown in figure 11 and compared with a calculated curve for 
6 i = 0.35° and with the small-deflection results of reference 3- The 
results calculated by the use of the large-deflection equations are in 
good agreement with experimental values and account for the frequency 
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increase measured in the vicinity of the buckling temperature. A sub- 
stantial improvement is noted over the predictions of small-deflection 
theory. In figure 12 the frequency ratio is shown as a function of 
temperature difference. 


CONCLUDING REMARKS 


An approximate analysis of the effects of initial twist and large 
deflections on the torsional stiffness of a cantilever plate subjected 
to nonuniform heating shows that for a perfectly flat plate the effective 
stiffness* and thus the torsional frequency* decreases with increasing 
thermal stress* just as predicted by small-deflection theory* and goes 
to zero when the buckling temperature difference is reached. Beyond 
the buc kli ng temperature* however, the stiffness* and thus the frequency, 
increases as the plate twists. If the plate has an initial twist* it 
begins to deform immediately upon heating and the stiffness decreases 
in much the same way as that of the perfect plate. The incremental 
stiffness of the initially twisted plate* however* reaches a minimum 
greater than zero before the theoretical buckling temperature difference 
is reached; further heating then increases the stiffness. The mi nl mum 
incremental stiffness is a function of the initial twist and, if the 
initial twist is sufficiently large* no reductions in stiffness are 
obtained and nonuniform heating then always increases the stiffness. 

These results for stiffness changes associated with the torsional 
frequency are also applicable to the stiffness of the plate with respect 
to small changes in the applied moment. In this case the applied moment 
has an effect similar to the effect of an Initial twist. If the applied 
moment exactly counteracts the Initial twist, the plate behaves in much 
the same way as the perfect plate. 

Calculated results were compared with available experimental data 
and were found to be in satisfactory agreement in view of the approximate 
nature of the calculations. 


Langley Aeronautical Laboratory* 

National Advisory Committee for Aeronautics* 
Langley Field* Va.* Ifey 15* 1957- 
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APPENDIX A 

DISCUSSION OF MODAL FUNCTIONS USED IN TEE ANALYSIS 


In the analysis of the torsional deformations of the cantilever 
plate by small-deflection theory, a deflection function containing six 
terms antisymmetrical in the coordinate y was. used. The function 
(eq. ( 8 ) in ref. 5 ) is 


v = An oX y + An l.X 


2J 


X X X 


The effect of various combinations of terms on the calculated critical 
temperature difference, torsional frequency, and moment-twist relation was 
investigated for a square cantilever plate (a/b = 2 ) by use of the small- 
deflection theory of reference 5* The following table shows the results 
obtained by starting with a single term A-^ and then progressing through 
other combinations to six terms : 


Terms used 

^cr 

Kl 

«o 2 

Ke 

M 

e 

A 12 

A i2 , A 22 

A 12-> a 22* A 32 

A 12> A 22 > A l4> ^2h 

A 12> A J2^ a 14* ^4' A j4 

103.9 
57. T 
50.4 
57-2 
50.0 

127.0 

77-5 

71.5 

77-4 

70.1 

0.324 x 10 6 
.278 

.273 

.278 

.273 


where and K 2 represent constants that include the plate dimensions 

and material properties. 

Inasmuch as the tabulated values were obtained from an application 
of the Rayleigh-Ritz procedure, the lowest value is the most accurate. 

The use of six terms improves the accuracy less than 2 percent over three 
terms and, consequently, is not worth the extra complication. The large- 
deflection analysis thus can be made by use of only three terms. The 
deflection function is nondimensionalized for this analysis and expressed 
in the form 

, + A ^(f) 3 + MSf] 

± 2 " 


1 + A-^ + A ^ 


(A2) 
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The relative values of these coefficients for the conditions given 
in the preceding table were obtained from, the small -deflection analysis 
of reference 5 for the square cantilever plate (a/b = 2) with the 
following results: 


Condition 

1 

A 1 

Ag 

1 + Ai + A 2 

1 + A]_ + A 2 

1 -f Ap + A 2 

Thermal buckling (A = 1, cd/cuq = 0) 

4.6l 

-5-75 

2.14 

Torsional, vibration (A = 0, c^/coq = l) 

3.63 

-4.06 

1.42 

Applied moment (A = 0) 

2.36 

- 2.01 

O .65 

Applied moment (A = O. 5 ) 

3.10 

- 3.20 

1.10 


The variation of fQ with x/a for these four conditions has been 
plotted in figure 13 . The table indicates a wide variation in the rela- 
tive values of the coefficients but the plotted data show that the deflec- 
tion shapes are all approximately the same. Then, so long as the calcu- 
lations are used to indicate changes from the initial conditions, any of 
the modes indicated should be satisfactory. 

In addition to the deflection function of equation (A2) , a stress 
function and a temperature function are required in the analysis and 
have been selected to correspond to those of reference 5- For convenience 
the functions are expressed in nondimens ional form as follows: 



f 3 - (*/ 

The exponent £ must be selected to describe best the measured tempera- 
ture distribution at the time of Interest. Note "that equation (A4) 
requires that the specified temperature difference AT be the difference 
between temperatures at the heated edges and the longitudinal center line. 
The function has been chosen so that the stresses vanish on the free 

boundaries. Along the root (x = 0) the function requires that the 
plate be free to expand parallel to the y-axis . Although this require- 
ment does not indicate a "built-in" condition, it probably resembles the 
actual test condition. 
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APPENDIX B 

EVALUATION OF INTEORAIS AND RELATED FUNCTIONS 


The Integrals indicated in the analysis have been evaluated in terms 
of the plate dimensions and the nondimensional functions given in 
appendix A. The integrals are as follows: 



12 



128 a t; 

15 (£ + 3) a + 1) 


(B2) 



2 



dx dy 




dx dy 


128 


15ah(l + Ai + A2) 


2\105 35 


+ £_ + 4- 1 

693 1 6 315 




20 / 






£ + 8A-, 



'd 2 f 2 3 2 f 2 | 

(dx 2 dy 2 


+ 2(1 



+ 3k Ag 2 + 8A 1 + + 24 a 1 a 2 + 

+ 6A x + + SA^ j 


dy dx 


(B4) 
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^5 = J J 5(x-a) {s(y-b) - 8(y-4b)jf 2 dy dx = 2 


(B 5 ) 


a 0 b 


*6 


■ a 


f 2 dy dx = 


gab / ■ 1 , . £l , g^2 . A l^ 

3(1 + A! + A^) 2 V 7 9 3 7 W 


(B6) 


For the square cantilever plate (a/b = 2) and the values of Aj_ 
and A2 given in appendix A for the thermal-buckling mode, the 
integrals I-j_ to I4 reduce to the following when Poisson's ratio p 
equals 0.33: 


t - f 128 ) 155 

1 W 2a 2 


(B 7 ) 


T . 256 i 

2 15 (£ + 3) (S + 1) 

(B8) 

1 

_ 1.520 

(B9) 


14 - 9 f 7 

(BIO) 



The integrals and Ig have not been included but are discussed 

subsequently. 

The critical temperature difference given by equation (l8) then 
becomes 


AT C r 


18.5 


D 

(5 + i)(S + 3) 

Bata 2 

5 


(Bll) 
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The twist and applied moment are related to their corresponding nondimen- 
sional quantities by equations ( 23 ) and (24) and may now he expressed as 


2 » 

cp = 0.035 Q\J^~ (B 12 ) 


m = 0.035 7 ^- =&- (B13) 

(GJ) 0 D 


Use of the thermal-buckling mode yields incorrect values for the 
natural frequency and the twist due to applied moment. In addition, the 
measured values differ slightly from the correct calculated values as a 
result of the imperfect clamping at the root of the cantilever plate. 

For these reasons, the measured initial values of torsional frequency 
and twist due to an applied moment are used in the calculations and the 
changes produced by the heating cycle are calculated from the following 
equations : 


A = 1 - cp^ (Bl4) 

= 1 - * + 3<P 2 “ (B15) 


The integrals 1^ and Ig are needed only to calculate the initial 

natural frequency and twist due to applied moment and, thus, have not been 
evaluated for the mode shape of thermal buckling. 

In the calculations, results of which are presented in figures 9 , 

10, 11, and 12, the following quantities were vised: 

E = 10.6 X 10 6 psi 
a = 12.8 X 10“ 6 °F” 1 

p = 0.35 
t = 0.25 in. 
a = 20 in. 
b = 10 in. 

8 i = 0.35° 
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Figure 1.- Dimensions and coordinate system of cantilever plate. 
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Figure 5-- Calculated incremental stiffness as a function of temperature 
ratio for several values of initial twist. 
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Pi 

Figure 6.- Calculated, ratio of minimum incremental stiffness to in itial 
incremental stiffness as a function of initial twist for several values 
of applied moment. 
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Figure 7.- Measured temperature histories 

line of the plate test 
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(a) M = 0 and 400 in -lb. 

Figure 10.- Comparison of measured and calculated values of twist as a 
function of temperature difference for three values of the applied 
moment . 
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(b) M = -400 in-lb. 


Figure 10.- Concluded. 




Figure 11.- Comparison of measured and calculated values of frequency ratio as a function of 

time for the rapidly heated plate. 
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Figure 12.- Comparison of measured and calculated values of the frequency 
ratio as a function of temperature difference. 









